N94-166?7 f 


A FIELD MEASURE OF THE "SHADE" FRACTION 


Alan R. Gillespie, Milton O. Smith and Donald E. Sabol 
Department of Geological Sciences, AJ-20 
University of Washington 
Seattle, WA 98195 


"Shade" has a technical definition peculiar to linear spectral mixture analysis of 
imaging spectrometer data: it is the reduction in radiance from a surface due to lighting 
conditions and geometry, and includes topographic shading described by photometric 
functions as well as shadowing at all scales. "Shade" is an important constituent of 
nearly all remotely sensed images, and is one endmember resolved in spectral mixture 
analysis, where it is represented as a fraction of the measured radiance and a characteristic 
spectrum. This spectrum is typically the null vector, provided the data have been cor- 
rected for atmospheric and instrument effects: i.e., "shade" is the radiance from an ideal 
black surface. 

In topographic shading, irradiance is reduced — typically in proportion to cos(i), 
where i (incidence angle) is the angle between the sun and the local surface normal vec- 
tors. Therefore, the radiance is lowered by a multiplicative factor. Shadowing occurs 
when i>90°, or when sunlight is blocked by adjacent high terrain; the only irradiance is 
down-welling skylight and bounce light from adjacent terrain. In spectral mixture anal- 
ysis, "shade" is regarded as an additive term. In this regard, it is an accurate description of 
the proportion of a scene that consists of ideal shadows ("checkerboard mixing"); how- 
ever, "shade" represents the multiplicative cos(i) factor as well, and here it should be 
interpreted as the proportion of shadow that would darken the scene an equivalent amount. 
In either case, the "shade" fraction is lessened by adjacency effects, because the scene has a 
non-zero reflectivity instead of the ideal black surface generally assumed. 

In spectral mixture analysis, field and laboratory reflectance spectra are utilized to 
represent endmembers other than "shade." Laboratory measurements are typically made at 
M)°, such that darkening due to shade is minimal. However, field measurements are 
typically made at greater incidence angles, and arc affected by both shading and shadows 
at a "subpixel" scale, due to roughness of the measured surface. This darkening leads to 
an underestimation of the scene reflectivity, and this underestimation is related to the 
"shade” fraction, F s , sought in spectral mixture analysis. Therefore it appears that field 
radiance data can be interpreted to yield a reflectance spectrum less affected by "shade," in 
greater agreement with laboratory measurements, and also to yield F s , which is a measure 
of scene roughness at the scale of measurement. 

Field spectra are commonly acquired by measuring the radiance from a small (10 1 
- 103 cm 2 ) area of natural surface, and comparing this to the radiance from a flat smooth 
standard, such as halon, of known reflectance and at the same viewing geometry. For the 
target, assumed to be Lambertian, 

R t (k) - S(X) cos(i) r t (k) + D (X) r (1) 

where Rt(V) is target radiance at wavelength S is the solar irradiance filtered by the 
atmosphere, r^ is the target reflectance, and D is the downwelling irradiance from the sky, 
assumed to be isotropic. When the halon or other standard is placed in front of the 
spectrometer, the measured radiance is 
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( 2 ) 


Rh(X) = S (k) cos(i) rh(X) + D(^) rh(X) 
where rj, is the halon reflectance. For S (X) » D(X.), it is assumed that 

r t (X.) = rh(?t) R t (X)/Rh(^)- 

However, in the visible spectrum D(X.) is variable but may be several percent of S(A,), so 
that this value of qCX) is in error. Furthermore, for a textured target "shade" must also be 
considered. Thus, 

R t (X) = S (k) cos(i) r t (^) (1 - F s ) + D(X) r^Jt). ( 4 ) 

It is clear that, even if S(K) » D(X), the apparent value of r t (X) is reduced in proportion to 
F s , such that 

q,(X) R t (X)/RhW - rt(X) (1 - F s ). (5) 


Equations (3) and (4) arc underdetermined, and it is not possible to calculate both 
r t and F s from them. However, it is possible to do this by making two additional 
radiance measurements, and Rfo\ in the field. For these measurements, the target and 
halon standard are both shadowed (for instance, by holding up a sheet of cardboard), such 
that S(A,) = 0: 

RhW = D(\) r h (X) ( 6 

R t W = D(X) r t (X). ( J 


Substituting Rh' and R t ‘ for Rh and R t in equation (3) produces a more correct 
estimate of r t . The "shade" fraction may be also calculated: 

F s = 1 - Rh’ (Rt - Ri) / (Rt’ (Rh - Rh')) ® 

where the X notation has been dropped for simplicity. 

It should be noted that the apparent value of F s will be independent of K 
provided the simplifying assumptions arc valid. An important refinement is to consider 
adjacency factors, which affect primarily the terms containing T>(X) but can have an 
impact of several percent on F$. We have incorporated a simplistic model of adjacency 
effects that assumes that: (1) the downwelling skylight is diminished according to i/n (i 
in radians), and (2) the measured radiance is increased by bounce light from adjacent 
terrain. It is assumed that the adjacent terrain is flat (i = 0) and partially shadowed 
according to some nominal value F^* of F s (c.g., 10%). The adjusted expression for 
target radiance is 

R t ~ r t S ( cos(i) (1 - F s ) + r t (1 - F s *) (i In) ) + r t D ( (1 - i/7i) + r t (1 - F s *) ). (9) 


We have calculated the sensitivity of F s to measurement error by evaluating the 
above equations. We have also investigated the effect of viewing geometry on F s and the 
magnitude of adjacency effects, using simulated surfaces having different roughness scales. 
This involved integration over a grid of 2500 cells, for which a digital terrain model was 
specified by a randum number generator. In addition, we have determined F s for 
constructed surfaces, using a CCD to measure radiance images. Finally, we have 
determined F s in natural field settings. 
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We discuss the use of F s in understanding VNIR images. An interesting 
ramification of this research is that F s is related to the surface texture at the subpixel 
level, a difficult characteristic to estimate otherwise. In vegetated terrains such as the 
Amazon rain forest, F s determined remotely may be an important parameter for 
estimating canopy architecture over large areas. It is also possible that the F s parameter 
may prove useful in relating VNIR and Radar images. 

Measurement of F s at a range of field scales (i.e., 10' 3 to 10° m) may provide a 
way to overlap with textural or roughness measurements made by microtopographic or 
stereometric surveys (10 -1 to 10 2 m). Such data have been used to characterize surface 
roughness as a fractal dimension for correlation with Radar backscatter coefficients and use 
in forward scattering models. 
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